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Abstract 

The expansion of an isolated hot spherical nucleus with excitation energy and its caloric curve 
are studied in a thermodynamic model with the SkM* force as the nuclear effective two-body inter¬ 
action. The calculated results are shown to compare well with the recent experimental data from 
energetic nuclear collisions. The fluctuations in temperature and density are also studied. They 
are seen to build up very rapidly beyond an excitation energy of ~ 9 MeV/u. Volume-conserving 
quadrupole deformation in addition to expansion indicates, however, nuclear disassembly above an 
excitation energy of ~ 4 MeV/u. 
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Understanding the density evolution of nuclear systems at moderate excitation energies 
(~ 2 — 8 MeV/u) is of much contemporary interest, both from the theoretical and the 


experimental point of view. This is particularly relevant ii^the context of nuclear multi- 

With excitation, the 


li LliC 


fragmentation at intermediate energy heavy ion collisions 
thermal pressure pushes the system towards expansion. At sufficiently high excitations, the 
system is ultimately driven towards the break-up density below which it ceases to exist in a 
mononuclear configuration and ultimately disassembles into many fragments. This has an 


Q 


of hot nuc 


ear systems; it is also 
in the astrophysical 


important bearing in modelling the equation of state 
of utmost importance in understanding explosive nucleosynthesis 
context. 

Experimental studies have been suggestive of decreasing break-up density with increasing 
excitation energy. Break-up densities have been determined from studies of correlation 


functions of emitted light particles from the source [3]. It is found that this density in general 
decreases with increasing excitation energy in the fragmenting system, however, the method 
of analysis leaves room for large uncertainties. QMD transport model calculations indicate 
that for medium-heavy systems, for excitation energy less than ~ 7 MeV/u, densities as low 


as ~ 0.35po can be reached 


where po is the ground state density. Break-up densities have 


also been determined from Coulomb barriers required to fit the intermediate mass ejectile 


spectra j9|, llO[ ; they have as well been determined from the analysis of apparent level density 
parameters required to fit the measured caloric curves At an excitation energy of ~ 8 
MeV/u, the deduced density comes down to ~ 0.3 — 0.4 po, but in the excitation energy 
range explored, the last two sets of results are not exactly in consonance; break-up densities 
derived from Coulomb barrier systematics are lower at higher excitations compared to those 
derived from caloric curve data. 


From a theoretical standpoint, statistical models [l|, |2, |3[ have been quite successful in 
explaining various observables related to nuclear fragmentation. The key parameter in these 
models is the freeze-out density p/. On reaching p/, the hot nucleus undergoes one-step 
prompt multifragmentation and the interaction among the generated fragments is assumed 
to be frozen out to change the fragmentation pattern. The freeze-out density is generally 
taken to be independent of the excitation energy. This density differs appreciably in different 


models. Whereas in the canonical or microcanonical models proposed in Refs. 


fl 


the 


density p/ is ~ 0.12 — 0.2 po, the corresponding density in the lattice-gas model 


liQ 
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is ~ 0.3 — 0.4 po- Such an uncertainty in both the theoretical and the experimental arena 
warrants a closer look at the relationship of the nuclear density with excitation energy. To 
explore the break-up density of a hot metastable mononuclear conhguration, Sobotka et al 


14| have recently performed a calculation where they utilise the fact that an isolated system 


with a given excitation energy pursues the maximal entropy conhguration for equilibrium. 
This calculation is schematic and involves a few ad hoc parameters; the results compare 


favorably with the data reported in Ref. 


Q. 


In the present communication, we address 


this problem with some basic microscopic inputs startin g w ith a Skyrme type effective two 


g w it 

m 


body interaction. We use the celebrated SkM* force [15ll, well known for its success in 
the description of a variety of ground-state properties of nuclei and diverse phenomena like 
nuclear hssion and nuclear collective modes. 

The experimental perspective in the laboratory is recalled briehy. When two nuclei collide 
at intermediate energy, a hot nuclear system (which may be initially somewhat compressed; 
the possible resulting collective how is ignored in this work) is formed which may be de¬ 
scribed statistically by an ehective temperature T. If the system were in a heat bath at the 
temperature T (canonical ensemble), the system could be described by equilibrium thermo¬ 
dynamics driving it to the minimum of the free energy. However, the system as prepared, 
is isolated with a hxed total excitation energy (microcanonical ensemble). The unbalanced 
thermal pressure induces expansion of the system in search of maximal entropy where the 
total pressure vanishes and the system is in equilibrium in a mononuclear conhguration. 
The energy of expansion is derived from the thermal energy, the temperature thereby de¬ 
creases. The density at this maximal entropy state is the lower limit (break-up density) 
for a mononuclear conhguration with a hxed excitation energy. The system may, however, 
gain further entropy from nuclear disassembly. A full-hedged study of the nuclear disassem¬ 
bly path is not taken up here, it is mimicked through a volume conserving deformation of 
the system at the various stages of expansion. For simplicity, we consider only quadrupole 
deformation. 

The initial state of the system is prepared by subjecting it in a heat bath at a chosen 


temperature T. Employing the SkM* interaction 


Q. 


this is done in a hnite temperature 


Thomas-Fermi framework within the subtraction scheme 0 10 . well suited for the de¬ 
scription of hot nuclei. The system so prepared is then detached from the heat bath and 
allowed to expand with constant total energy in pursuit of the maximal entropy state. The 
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expansion is simulated through a self-similar scaling approximation for the density: 

px{r) = AV(Ar), (1) 

where the scaling parameter A lies in the range 0 < A < 1 and p{r) is the base density prohle 
generated in the Thomas-Fermi procedure. 

In actual calculations, we hrst £x the total excitation energy E*. The base density prohle 
of the system is generated at a chosen temperature T such that the excitation energy for this 
density prohle is less than the given E*. The system is then allowed to undergo a self-similar 
expansion till the total excitation energy (thermal plus expansion) reaches E* at some value 
of A < 1. For any density prohle, the excitation energy is calculated as 


E* = E{X,T) - E{X = 1,T = 0) (2) 

from the Skyrme energy density functional The corresponding entropy for the ex¬ 

panded conhguration S{X,T) is computed. The calculations are repeated for diherent T; 
the conhguration corresponding to the maximum of the entropy prohle S{X, T) is the desired 
equilibrium conhguration at the energy E*. 

In the subtraction procedure nL Il 7l| the base density for the hot nucleus is given by p(r) = 
Pngij) — Pgij)) where Pg{r) is the density of the surrounding gas representing evaporated 
nucleons in which the system is immersed to maintain equilibrium at the temperature T and 
Pngij) is the density of the nucleus-plus-gas system. The density prohle p{r) of the heated 
nucleus is then independent of the size of the box in which calculations are done, the density 
and pressure being zero at large distances. In the Thomas-Fermi method, the densities pi{r) 
for neutrons or protons [i stands for ng or g) are given by 


PiPr) = 


27r2 


2m* (r) 


'Vi{r) 


x/e- Vi{r)f{e, p, T) de. 


( 3 ) 


Here m*{r) is the ehective fc-mass of the nucleon, Vi{r) is the single-particle potential, / is 
the Fermi occupation factor, and p is the chemical potential which is same in both the ng and 
g phases. The ehective mass, single-particle potential, and the chemical potential are isospin 
dependent. The chemical potentials are determined from particle number conservation (iV 
is the neutron or proton number): 


N= / g{e,T)f{e,p,T)de, 


( 4 ) 
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where the sing 
is given by 


e-particle level density g{e, T) of the hot nucleus in the subtraction procedure 
4^2 


m 

9{e,T) = 


71 


(mng)^ \Je-Vng{r) - (m^^ sJe-Vg{r) 


r dr. 


(5) 


The effective mass, single-particle potential, and the chemical potential are evaluated at the 
appropriate temperature and scaled densities, so also the single-particle level density and 
the entropy. The latter reads as 


S{\T) = - / 9(£, D 1/ 111 / + (l-/)lii(l-/)]&. 


( 6 ) 


The total entropy is the sum of the neutron and proton contributions. 

For our study, we have chosen ^®°Sm as a representative system. The central density pc of 
the calculated mononuclear equilibrium conhguration in units of po (the central density of the 
unexpanded nucleus at T = 0) |l^ is displayed in Fig. 1 as a function of the excitation energy 
and compared with the experimental data. The thin (dashed and solid) lines correspond 
to the canonical ensemble calculations, he., when the system is in a heat bath; the thick 
lines are those for the isolated expanded nucleus at equilibrium. The hlled circles are the 
experimental points extracted from the apparent level density parameters for the mass 
selection 140 < A < 180, where A is the mass number of the system, and the empty squares 


are the ones obtained from Coulomb barrier systematics 


10| for Au-like systems. 


The effective mass m* dehned previously comes from the momentum dependence of the 
single-particle potential, which is the fc-mass rrik- However, m* should have a frequency 
dependent mass-factor m^/m: 

TTi/f. , , 

m = m -. (7) 

m m 

The a;-mass originates from the coupling of the single-particle motion with the collective 
degrees of freedom. This has the effect of bringing down the excited states from high energy 
to lower energy near the Fermi surface, thus increasing the many-body density of states 
[m^/m > 1, see Eq. (8)] at low excitations. It may a priori have a signihcant role to play in 
the present context as the system can accommodate comparatively more entropy at a given 
excitation energy. The self-consistent evaluation of the ci;-mass is beyond the scope of the 
present work; we use the phenomenological form jl^. 1^ 


m 


= 1 — 0.4 As exp 


T 


21A-3 


1 dpipr) 
p(0) dr ’ 


( 8 ) 


5 
















where p{r) is the density prohle at temperature T. The effect of is incorporated 00 
by replacing the single-particle potential V with {m/mS)V in Eqs. (3)-(5). 

In Fig. 1, the solid (dashed) lines refer to calculations performed with (without) inclusion 
of m^/m. With increasing temperature, m^/m tends to unity. We find that the values of pc 
calculated at moderate or higher excitations with and without the inclusion of the a;-mass are 
practically the same; however, the relative densities pdpo for the two cases run parallel there 
because the values of po are different in the two situations. The calculations for the canonical 
ensemble (thin lines) terminate at an excitation energy of ~ 5.5 MeV/u corresponding to 


T ~ 8.5 MeV, the limiting temperature 


00 


for this nucleus beyond which it is unstable 


in a heat bath. The thick curves, as compared to the thin ones, show appreciably lower 
values of density that compare reasonably with the densities extracted from the analysis 
of the caloric curve measurements performed in 11| (filled circles). The importance of a 
microcanonical treatment is thereby indicated. The fit with the data obtained from Coulomb 
barrier systematics (filled squares) 0| is relatively poor; however, an ambiguity in their 
extraction procedure has been pointed out recently |22| . 

The (u-mass does not appear to have a very distinctive role to the density evolution with 
excitation energy. In the excitation range indicated in the figure, it turns out that with 
inclusion of the cn-mass, the equilibrium configuration corresponds to a lower value of the 
temperature as compared to that obtained with this tends to increase the central 

density pc in the former case. On the other hand, the scale parameter A is found to be 
comparatively lower with inclusion of the a;-mass, resulting in a lower value of pc- The 
combined effect of these opposing tendencies then results in the near equality of pc in both 
the calculations. 

The expansion of the nucleus has an important bearing on the correlation of the excitation 
energy with temperature. The caloric curve so obtained for the expanded hot nucleus 
in equilibrium is displayed in the upper panel of Fig. 2. At a fixed excitation energy, the 
system cools down with expansion and therefore the recorded temperature at the equilibrium 
configuration is significantly lower than that for the unexpanded nucleus prepared initially 
with the same excitation. In the figure, T refers to the canonical temperature. We have 
checked that the microcanonical temperature obtained from T~^ = dSe^/dE*, where S'eq is 
the total entropy and E* the excitation energy at equilibrium, is not much different from 
the canonical one. For comparison, a representative set of experimental data for medium 
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mass nuclei jSj are also shown in the hgure. The notations used have the same meaning 
as in Fig. 1. The effect of isolation is not appreciable at relatively lower excitations, but is 
more manifest at higher excitations. The calculated microcanonical results are seen to agree 
nicely with the experimental caloric curve. The canonical caloric curves with and without 
the inclusion of the a;-mass merge at high temperatures because rrii^/m does not differ much 
from unity. The microcanonical caloric curves, on the other hand, even at high excitations 
do not merge; the highest temperature encountered in these calculations is ~ 8 MeV where 
rrii^/m is typically ~ 1.02 at the surface region affecting the equilibrium configuration a little 
which is amplified in the caloric curve. 

At a given temperature, the system is found to have more excitation with inclusion of the 
a;-mass as more states are available near the Fermi surface to absorb more energy. Beyond 
E*/A ~ 6 MeV, T seems to saturate and at ~ 8 MeV, a downward slope in the caloric 

curve is apparent for the isolated system implying negative heat capacity. The expansion 
energy [F^expn = E{\, T) —E{X = 1, T)] comprises a significant part of the total excitation at 
higher values of E* inducing the above-mentioned characteristics in the caloric curve. The 
expansion energy is found to be in good agreement with that obtained in the Expanding 


Emitting Source model of Friedman as shown in the lower panel of Fig. 2 The filled 
diamond refers to an experimental estimate [^, which is close to the prediction from our 
calculations. 

The entropy profile obtained at a fixed excitation energy may be employed to calculate the 
mean values and the fluctuations around the mean of the temperature and of the density or 
volume. The probability of finding a configuration with scale parameter A and temperature 
T at a fixed excitation energy E*/A is given by 

W{X,T) ^ (9) 


where S{X,T) is the total entropy (jHI) of the given configuration. The n-th moment of the 
central density pc is then given by 


(p:(A,T)) 


J e^(^E)pn^X,T)dp 

JeS(^E)dp ’ 


( 10 ) 


which allows to calculate the mean {pc) and the variance cr^ = {pi) — {pc)"^- Similarly, the 
mean and variance of the temperature at a fixed excitation can be evaluated. 

For a thermodynamic system, the average and the most probable (equilibrium) value of an 
observable are the same. For a finite system, however, they may be different. Experimentally, 
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the average value is the relevant quantity. The average and the most probable values of 
the temperature and the specihc volume Vc (= 1/Pc) in units of vq (= 1/po) for the system 
considered are displayed in the upper and lower panels of Fig. 3. It is seen that the differences 
between the averages and the most probable values of the temperature and the break-up 
volume (or density) are not very signihcant. The fluctuations in temperature and the 
specihc volume (measured in units of vq) are shown in Fig. 4. The huctuations rise smoothly 
up to the excitation energy E*/A ~ 9 MeV, beyond which this build-up is very sudden; this 
is particularly more pronounced for the volume huctuations. This large density huctuation 
indicates that beyond E*/A 9 MeV the system becomes unstable and breaks up in many 
pieces. It turns out that the negative branch of the specihc heat (Fig. 3) and the large 
huctuations start at around the same excitation energy. A possible close correlation between 
them is thereby indicated. 

So far, for the search of the maximum entropy conhguration, the shape of the excited 
expanding system has been constrained to a spherical one; a possible deformation path 
along with expansion might also contribute additional entropy and mimics a fragmentation 
channel. To investigate this aspect, at all stages of the expansion at a hxed excitation 
energy, a volume conserving quadrupole deformation is explored. In a volume conserving 
deformation, only the surface and Coulomb energies change. To calculate these changes, a 
sharp surface approximation to the density prohle is made (i?sharp = which also 

facilitates the calculation of entropy from deformation. At a deformation /?, the excess 


Coulomb energy of t 
/(/9) is given in Ref. 


le expanded deformed nucleus is SEc{X, (3) = F^c(A, 0)/(/3); the function 
The surface free energy due to deformation is given by SE{X, j3) = 
5M(A, /3)cr(/9, T), where SA is the excess surface area of the nucleus arising out of deformation. 

y p and temperature T is taken as a{p,T) = 


The surface tension coefficient a at a densi 
a(T)g{p). The temperature dependence [2^ of the surface tension is given by a(T). The 
surface tension has its maximum value at the ground-state density; for an expanded system, 
a decreases which can be well represented by g{p) taken to be a polynomial in p. We 


ci Q 


have calculated this density dependence from the prescription of Myers and Swiatecki 
using the scaling approximation to the ground-state density prohle of semi-inhnite nuclear 
matter. The excess entropy from deformation is then calculated as 5S = —d{5E)/dT\p 
which immediately gives the excess surface energy 5F^surf = + T5S. The total excitation 

energy of the expanded deformed system is then T*(A, {3) = £*{X, 0) -|-(5i?c(A, (3)+6Esnri{X, (3) 






and the corresponding entropy is S'(A,/3) = S{X, 0) + 6S{X,j3). 

Along the deformation path, a barrier is faced coming from the interplay of the Conlomb 
and snrface energies; it decreases with increasing temperatnre and increasing expansion. In 
actnal calcnlations, the scale parameter A is adjnsted snch that for a chosen temperatnre, 
the excitation energy at the top of the barrier S*{X,f3) matches the given excitation E*/A. 
This is repeated for different temperatnres and the maximnm entropy among these different 
conhgnrations is selected. If this entropy exceeds that for the expanded spherical eqnilibrinm 
conhgnration, then the deformed shape is favonred leading to the fragmentation channel. 

This extra entropy gained dne to deformation (we call this AS') over that at the spherical 
eqnilibrinm shape is shown in the top panel of Fig. 5 as a fnnction of excitation energy. 
For E*/A less than ~ 4 MeV, AS is negative in onr chosen restricted deformation space 
and so a spherical eqnilibrinm conhgnration is more probable. Above this excitation energy, 
fragmentation resnlting from deformation is more favonrable. Beyond E*/A ~ 9 MeV, the 
barrier vanishes and the system nndergoes spontaneons fragmentation from the spherical 
eqnilibrinm conhgnration. 

In the middle panel of Fig. 5, the eqnilibrinm central densities with (solid line) and with- 
ont (dashed line) deformation are displayed as a fnnction of E*/A. With deformation, the 
maximal entropy conhgnration occnrs at a lower temperatnre with a smaller scale param¬ 
eter A resnlting in a somewhat rednced density as seen. The corresponding caloric cnrves 
are shown in the bottom panel. All the above calcnlations pertain to prolate deformation, 
oblate shapes are fonnd to have lesser entropy. 

We have addressed to some gross featnres of a mononnclear conhgnration at medinm and 
high excitations in a semi-microscopic framework based npon a realistic ehective nnclear 
interaction. The calcnlated break-np densities for the nnclens in the micro canonical for- 
mnlation are i n g eneral in good agreement with the ones extracted from the experimental 


re inge: 

hs [111. 


data analysis [ll|. The generated mononnclear caloric cnrve also compares very well with 
the experimental results. The plateau observed in experimental data has been taken to 
be suggestive of a possible phase coexistence j22|. In our calculations, the plateau comes 
naturally from nuclear expansion with excitation. For the bloated spherical mononuclear 
conhgnration, the rapid build-up of huctuations, particularly in the density, is suggestive of 
the instability of this conhgnration against prompt multifragmentation beyond E* jA ~ 9 
MeV. Expansion with deformation degrees of freedom may have signihcant ehects on the 
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physical observables; to have an orientation on the role of deformation, we have considered 
volume-conserving quadrupole deformation. It is found that above E*/A ~ 4 MeV the sys¬ 
tem favours deformation, a precursor to fragmentation. The global trends provided by the 
present model, particularly with the inclusion of the frequency dependence in the effective 
mass, are qualitatively consistent with the experimental data as well as with the 


ones obtained from other mononuclear formulations |l4128| . At high excitations, the collec¬ 
tive flow may play a signihcant role in the nuclear collision process. This has been ignored 
in the present calculation; it would be worthwhile to investigate its effect on the break-up 
density and also on the caloric curve. 
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Figure Captions 


Fig. 1 The calculated equilibrium density in units of po as a function of the excitation energy 
per nucleon is compared with the experimental data. For notations, see text. 

Fig. 2 The mononuclear caloric curve (upper panel) and the expansion energy per nucleon 
(lower panel). 

Fig. 3 The average and the equilibrium values of temperature (upper panel) and the specihc 
volume (lower panel) of the nucleus as a function of excitation energy. 

Fig. 4 The variance in the temperature (upper panel) and the specihc volume (lower panel) 
as a function of excitation energy. 

Fig. 5 The deformation entropy AS (top panel), the equilibrium density in units of po (middle 
panel) and the caloric curve (bottom panel). The solid lines represent the results 
with deformation and the dashed lines correspond to the spherical conhguration. The 
experimental points (hlled circles and open squares) are the same as those given in 
Figs. 1 and 2. 
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heat bath, m^/m=1 

- heat bath, Eq.(8) 

-Isolated, mym=1 

- Isolated, Eq.(8) 

• Expt. [11] 

□ Expt. [9,10] 
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